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We report a model-independent measurement of the entropy, energy, and crit- 
ical temperature of a degenerate, strongly interacting Fermi gas of atoms. The 
total energy is determined from the mean square cloud size in the strongly 
interacting regime, where the gas exhibits universal behavior. The entropy 
is measured by sweeping a bias magnetic field to adiabatically tune the gas 
from the strongly interacting regime to a weakly interacting regime, where 
the entropy is known from the cloud size after the sweep. The dependence of 
the entropy on the total energy quantitatively tests predictions of the finite- 
temperature thermodynamics. 

Strongly interacting Fermi gases are of great interest, as they exhibit universal thermody- 
namic behavior, where the properties are independent of the details of the microscopic interac- 
tions ([7] [2] El®. These gases provide models for testing nonperturbative many-body theories in 
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a variety of fields from neutron stars and nuclear matter (|5]|6][2][7]) to quark-gluon plasmas © 
and high temperature superconductors ©. Hence, thermodynamic experiments on strongly 
interacting Fermi gases are of great importance. 

In studies of the thermodynamics of these systems, where thermometry is difficult (1701) . 
entropy measurement plays a central and fundamental role. We report the measurement of the 
entropy S of a strongly interacting Fermi gas as a function of its total energy E. The results yield 
the temperature T via the elementary thermodynamic relation 1/T = dS/ dE. Our experiments 
quantitatively test recent predictions of the entropy based on microscopic many-body theory, 
yield the dependence of the energy on temperature, and determine the critical temperature for 
the superfluid transition without invoking any specific theoretical model. 

Strongly-attractive Fermi gases exhibit both fermionic and bosonic features, and have been 
studied intensely for several years in theory (Iiilli2ll73l l9l) and experiment (I71 l74lli5lli6ll77lli#lli9l) . 
Measurements of the heat capacity (l20l) and collective mode damping versus energy (|27|) reveal 
transitions in behavior, which have been interpreted as a superfluid transition in this system (1201) . 
Recently, the observation of vortices ([22]) has provided a definitive proof of a superfluid phase. 
However, there have been no model-independent studies of the thermodynamic properties. 

A strongly interacting Fermi gas is prepared using a 50:50 mixture of the two lowest hy- 
perfine states of 6 Li atoms in an ultrastable CO2 laser trap with a bias magnetic field of 840 
G just above a broad Feshbach resonance at B = 834 G (|23l) . The gas is cooled to quantum 
degeneracy by lowering the trap depth by a factor of ~ 1000 ©. Following forced evaporation, 
the trap depth Uq is recompressed toUo/ks = 10 /iK, which is large compared to the energy per 
particle of the gas. Here ks is the Boltzmann constant. After this procedure, the initial energy 
is close to that of the ground state, as described below. 

At the final trap depth, the measured trap oscillation frequencies in the transverse directions 
are uj x = 2n x 670 Hz and uj y — 2tt x 760 Hz, while the axial frequency is cu z = 2ir x 30 
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Hz @ 840 G and lu z = 2ir x 32 Hz @ 1200 G. Note that axial frequencies differ due to 
the small change in the trapping potential arising from the bias magnetic field curvature. The 
total number of atoms N ~ 1.3(0.2) x 10 5 is obtained from absorption images of the cloud 
using a two-level optical transition at 840 G. The corresponding Fermi energy E F and Fermi 
temperature T F for an ideal (noninteracting) harmonically trapped gas at the trap center are 
E F = k B T F = hujiy'iN) 1 ^ , where uj = (co^^cll.) 1 / 3 . For our trap conditions, we obtain 
Tp ~ 1.0 ^K. 

The total energy per particle, E, of the strongly interacting gas is measured in a model- 
independent way from the mean square size in the axial direction (4). In this strongly interact- 
ing regime, the zero energy s-wave scattering length a$ is large compared to the interparticle 
spacing, which is large compared to the range of the two-body interaction, so that the gas is uni- 
versal (|2][7]|2]). Then, the local pressure is P = 28/3, where 8 is the local energy density (ji]®. 
Using force balance for a trapping potential U, VP + riWU = 0, where n is the local density, 
one then obtains the total energy per particle E = 3mw 2 z (z 2 )s4o (1 — «) or 

|-=^(l-K), (1) 

t,p z F 

where (-2 2 ) 8 4o is the mean square axial cloud size measured at 840 G and m is the 6 Li mass. 
Here, zp is defined by ?>muj\ zp = Ep, and is weakly dependent on the magnetic field through 
the trap frequencies. The correction factor 1 — k arises from anharmonicity (|2?)) in the shallow 
trapping potential U ^ 10 E F used in the experiments. We find that k varies from 3% at our 
lowest energies to 13% at the highest. 

The entropy of the strongly interacting gas at 840 G is determined using an adiabatic sweep 
of the magnetic field to a relatively weakly-interacting regime at 1200 G, where a reference 
entropy can be estimated from the mean square axial cloud size (,2 2 )i2oo- At 1200 G, as = 
-2900 bohr (23), and k F a s = -0.75 for our shallow trap, with k F = ^2mk B T F /h 2 . At 
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kpa = —0.75, we expect that the dependence of the entropy on the cloud size should be close 
to that of an ideal noninteracting Fermi gas with primarily a small mean-field reduction in 
the ground state cloud size. This conjecture is supported by the observed ballistic expansion 
of the cloud at 1200 G, even at our lowest temperatures, which shows that the gas is nearly 
normal. We also find that the calculated ideal gas entropy differs from a many-body result 
for kpCbs = —0.75 (|25l) by less than 1% over the range of energies we studied, except at the 
point of the lowest energy, where they differ by 10%. For this comparison, we slightly shift the 
ground state size of the ideal gas to coincide with that calculated for kpCLs = —0.75. Hence, the 
reference entropy at 1200 G is nearly identical in shape to that for an ideal gas, and provides a 
model-independent estimate of the entropy of the strongly interacting gas. 

Ideally, a sweep from 840 G to a magnetic field of 528 G, where the scattering length van- 
ishes, would produce a noninteracting gas (k F a$ = 0), where the entropy is precisely known. 
Unfortunately, adiabatic formation of molecules (26) and subsequent molecular decay at fields 
below resonance (17) cause unwanted heating for such a downward sweep. 

To measure the entropy as a function of energy, we start with an energy near the ground 
state and controllably increase the energy of the gas by releasing the cloud for an adjustable 
time and then recapturing it, as described previously (20). After recapture, the gas is allowed to 
reach equilibrium for 0.7 s. This thermalization time is omitted for measurement of the ground 
state size, where no energy is added. 

After equilibrium is established, the magnetic field is either ramped to 1200 G over a period 
of 1 s, or the gas is held at 840 G for 1 s. In either case, after 1 s, the gas is released from 
the trap for a short time to increase the transverse dimension of the cloud for imaging, without 
significantly changing (less than 0.5%) the measured axial cloud size. 

We find that the magnetic field sweep is nearly adiabatic, since the mean square size of 
the cloud at 840 G after a round-trip- sweep of 2 s duration is found to be within 3% of that 
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obtained after a hold time of 2 s at 840 G. However, we also find for our shallow trap that there 
is a magnetic field and energy independent heating rate, which causes the mean square size 
to slowly increase at a rate of (z 2 ) = 0.024 zp/s, corresponding to 24 nK/s in energy units. 
Since we desire the energy and entropy just after equilibration, we subtract (z 2 ) x 1 s from the 
measured mean square axial dimensions for both the 840 G and 1200 G data. The maximum 
correction is 5% at the lowest energies. 

Fig. Q] shows the ratio of the mean square axial cloud size at 1200 G (measured after the 
sweep) to that at 840 G (measured prior to the sweep), as a function of the energy of the strongly 
interacting gas at 840 G. The energy at 840 G is directly measured from the axial cloud size at 
840 G using Eq. [Q The displayed ratio and energy scale are independent of the atom number 
and trap parameters. This is accomplished by measuring the mean square sizes at each field in 
units of z F for the given field and atom number. The total data comprise 900 measurements 
which have been averaged in energy bins of width AE = 0.04 Ep. 

The red solid line shows the predictions obtained by equating the entropies calculated at 
1200 G and near resonance (25). The predicted curve exhibits a rapid drop followed by a 
slower decline to unity, in very good agreement with the data in the low and high energy re- 
gions. However, the data deviate significantly from the prediction in the region centered near 
E — E ~ 0.4 Ep, where the entropy changes behavior as described below. 

We note that potential energy has been measured previously in 40 K (|27l) at a Feshbach res- 
onance and after an adiabatic sweep to the noninteracting regime. In Ref. (|27l) . the resulting 
potential energy ratios are given as a function of the temperature of the noninteracting gas. 
In contrast, by exploiting universality, our cloud size ratios are referred to the total energy in 
the strongly interacting regime, which enables a measurement of S(E) and T for the strongly 
interacting gas. 

For our measurements of S(E), the origin for S = is determined by the cloud sizes (z 2 ) 
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for the ground states at 840 G and 1200 G. These sizes are estimated from the data at the lowest 
temperatures. In the harmonic approximation, the ground state obeys (z 2 )o/zp = (3/4) 
where £ = 1 + [3 is the ratio of the energy per particle of the strongly interacting gas to that 
of a noninteracting gas with the same density (|5]|6]|2][7]). At our lowest temperatures, including 
anharmonicity arising from the gaussian trapping potential U tr a P , we find /3 e // = —0.50(0.04) 
at 840 G. For our trap parameters, this corresponds to (3 = —0.54(0.04) at 834 G, using the 
estimate of Ref. Our result is in good agreement with recent measurements based on 

the axial cloud size, where (3 = -0.54(0.02) (HB, (3 = -0.54(+0.05/ - 0.12) (H3) and with 
recent calculations, (3 = -0.56 ©, (3 = -0.545 CM)), (3 = -0.564 $28$. Using our measured 
(3 e ff = —0.50, the ground state energy per particle for the strongly interacting gas is (|20|) 
E = (3/4)y/£E F , yielding E = 0.53 E F and (z 2 ) /z 2 F = 0.55 at 840 G. 

We can predict the ground state cloud size at 1200 G using the equation of state at zero 
temperature. An approximate equation of state for the chemical potential versus local density, 
fj,(n), is given in Ref. (|28l) . Very good agreement with quantum Monte Carlo calculations is 
obtained for negative scattering lengths, which is the region of interest to us. We invert the 
equation of state to find n(/i), and then using /i = fi g — U tra p, we determine the density for 
a gaussian potential U tTap to include anharmonicity. Normalization to the number of atoms 
yields the global chemical potential fx g and the mean square cloud size. At 1200 G, where 
kpas = —0.75, we find (z 2 )q/z f = 0.69. Our measurements at the lowest temperatures yield 
(z 2 )q/z 2 f = 0.72(0.02) at 1200 G, in agreement with the calculated value. Hence, at both 1200 
G and 840 G, we obtain clouds nearly in the ground state and the corresponding cloud size ratio 
0.72/0.55 = 1.31 shown in Fig. [B 

To convert the data of Fig. Q] into an entropy measurement, we calculate the entropy at 1200 
G as a function of the ratio ( (z 2 ) — (z 2 )o)/ zp, which is determined from the axial cloud size data 
at 1200 G. This method automatically assures that S = corresponds to the measured ground 
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state (z 2 )o at 1200 G, and compensates for small shifts between the calculated and measured 
ground state sizes. Then, S[((z 2 ) — (z 2 ) )/z F ) is obtained from a many-body calculation at 
kpds = —0.75, assuming an isotropic gaussian trapping potential, which automatically corrects 
for anharmonicity (|25ll37T) . As discussed above, nearly identical results are obtained if we 
assume that the entropy at 1200 G is that of an ideal Fermi gas in the same potential. 

Fig.[2]shows the entropy (blue dots) of the strongly interacting gas at 840 G as a function of 
its energy in the range < (E — E )/Ep < 1.4. The maximum energy is restricted to avoid 
evaporation in the shallow trap, which can reduce the energy and the atom number during the 
time of the magnetic field sweep. The entropy of the strongly interacting gas differs significantly 
from that of an ideal gas (lower orange dot-dash line), which has a larger ground state energy 
E I0 = 0.75 E F , so that E I0 — E = 0.22 E F . To compare the curve shape for the measured 
entropy to that of an ideal gas, the ideal gas entropy is also plotted with its energy origin shifted, 
so that S = at E — E = (upper orange dot-dashed line). In addition, the data are compared 
to predictions in the resonant regime based on pseudogap theory (|25][37]) (dotted red line) and 
quantum Monte Carlo methods (dashed green line) (|32ll33l) . 

The temperature is determined in a model-independent manner from l/T — OS/dE. This 
requires parameterizing the S(E) data to obtain a smooth curve. The simplest assumption 
consistent with S(E = E ) = is to approximate the data by a power law in E — E . However, 
one expects that below and above the superfluid transition at a critical energy E c , the power law 
exponents will be different. This suggests the simple form, 



where the fit parameters are a, b, d and E c . A fit with this parametrization yields a x 2 P er degree 
of freedom ~ 1, a factor of 2 smaller than that obtained by fitting a single power law to all of the 
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for < E - E < E c 




c; 



(2) 



data. However, Eq. [2] ignores the smooth transition in slope near E c , as required for continuity 
of the temperature, since the detailed critical behavior near E c is not resolvable in our data. 

Fitting the data of Fig. [2] with Eq. [2l the critical energy is found to be (E c — E )/E F = 
0.41 ± 0.05, with a corresponding critical entropy per particle S c = 2.7(±0.2) k B . Below E c , 
the entropy varies with energy as S<(E) = k B (4.6 ± 0.2) [(E - E ) /E F } om±om . Above E c , 
we obtain S > (E) = k B (4.0 ± 0.2) [(E - E ) /E F ] OA5±om . We find that the variances of a and 
b have a positive correlation, so that S(E) is determined more precisely than the independent 
variation of a and b would imply. The change in behavior near E c is shown clearly in the inset 
of Fig. [2] and in the log-log plot of Fig. [3] 

The power law exponent below E c , b = 0.61, falls between that of an ideal harmonically- 
trapped Fermi gas, where a Sommerfeld expansion at low energy yields S oc (E — Eq) 1 / 2 and 
that of an ideal harmonically-trapped Bose-Einstein condensate, where S oc (E — E ) 3 ^ 4 . By 
contrast, above E c , the exponent d = 0.45, is close to the result we obtain by fitting a power 
law to the entropy of an ideal gas, i.e., Sj(E — E I0 ) oc (E — E I0 ) b . In this case, b = 0.485 
for E — Ej below 0.41 Ep and b = 0.452 above. This is consistent with the cloud size ratios 
shown in Fig. CD which converge to unity at higher energies. 

The fit parameters from the data can be compared to those obtained from fits to the theo- 
retical curves shown in Fig. [2l The pseudogap theory (|25]|37l) predicts E c — E = 0.36 Ep, 
and S c = 2.16 Nk B . Using Eq. [2] to fit the theoretical curve below the predicted E c , we find 
S<(E) = k B (4.244 ± 0.003) [(E - E )/ E F ] omi±om5 . For the quantum Monte Carlo treat- 
ment <|22]|22J, whicn P redicts Ec - E = 0.32 E F , we find S c = 2.17 Nk B , and S<(£) = 
k B (4.35 ± 0.05) [(E - E ) /E F ] 0S13±0m7 . Here the error estimates do not include the error in 
the theoretical curves. The small variances indicate that the power law fit closely approximates 
the theory, showing that Eq.[2]is a reasonable parametrization. 

The energy versus temperature E(T) is determined from the derivative of the fit function 
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S(E). For E < E, 



E — Ei 



E F 



■o 




(3) 



From the best fit to the entropy data, where a = 4.6, b = 0.61, E c = 0.41 E F , we obtain 



We estimate the critical temperature T c using the measured value of E c — E = (0.41 ± 
0.05) Ep. Here, we interpret E c as the critical energy for the superfluid transition. We note 
that using E = 0.53 Ep yields E c = (0.94 ± 0.05) Ep. This value is consistent with our 
previous measurements based on the heat capacity, where we observe a change in behavior at 
E = 0.85 E F (20), and in collective mode damping (1271) . where a plot of the damping rate versus 
energy (rather than empirical temperature) shows a change in behavior near E = 1.01 E F . 

Ideally, to obtain T c , the fit S(E) should have a continuous slope near E c . Since our fit 
function has different slopes above and below E c , we approximate the true slope by the average, 
as expected for the tangent to a smooth curve. Inverting Eq. [3] yields T/T F = 0.36 [(E — 
E )/E F } - 39 and T C< /T F = 0.25. Similarly, for E(T) > E c , we find T/T F = 0.56 [(E - 
E )/E F } - 55 and T C> /T F = 0.34. Assuming that 2/T c ~ 1/T C< + 1/T C> , we find T C /T F = 
0.29(0.02). Here, the error estimate includes the cross correlations in the variances of a, b, E c , 



The measured critical temperature T C /T F = 0.29(0.02) can be compared to our previous 
estimate of T C /T F = 0.27 from an experiment with a model dependent temperature calibra- 
tion (20). Moreover, the result 0.29 is in good agreement with predictions for trapped atoms, 
0.29 <H3>, 0.30 423), 0.31 (123), 0.30 d221), 0.26 C2J and 0.27 

Transition temperatures also have been predicted for a uniform gas, ksTc/ep = 0.152 (|33|) 
and kBT c /e* F = 0.160 (l?7l) . These also can be compared to our measured T c . Here is the 
Fermi energy corresponding to the uniform density. By contrast, we determine the ratio T C /T F , 
where T F is the Fermi temperature for a noninteracting gas at the center of a harmonic trap. If 



(E — Eq)/ E F = 14 (T/Tp) 



2.56 



and d. 
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we assume that e F corresponds to the central density of the strongly interacting gas in our trap, 
then £e* F = ^k B T F <|20J. From this, we estimate T C /T F = k B T c /(e* F ^). For Ref. 
we assume £ = 0.44 (7), and obtain T C /T F = 0.23. Ref. © calculates f = 0.36 yielding 
TJT F = 0.27. 
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Figure 1: Ratio of the mean square cloud size at 1200 G, (z 2 )i2oo to that at 840 G, (z 2 )s4o- The 
data is obtained by adiabatically sweeping a bias magnetic field from 840 G, where the Fermi 
gas is strongly interacting, to 1200 G where it is weakly interacting. E 8i0 is the total energy of 
the strongly interacting gas at 840 G prior to the sweep, E is the ground state energy at 840 G, 
and E F the Fermi energy of a noninteracting gas. The solid line shows the theoretical prediction 
based on the calculated entropies (l25i The ratio converges to unity at high energy, as expected 
(dashed green horizontal line). 
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Figure 2: Measured entropy of a strongly interacting Fermi gas at 840 G versus its total energy 
(blue dots). The entropy is estimated from the measured cloud size at 1200 G after an adiabatic 
sweep of the magnetic field from 840 G. Lower orange dot-dashed curve- ideal gas entropy; 
Upper orange dot-dashed curve- ideal gas entropy with the ground state energy shifted to E ; 
Red dots- pseudogap theory (25); Green dashes- quantum Monte Carlo prediction (132)) . Inset- 
entropy versus energy data showing knee at E c — E = 0.41 E F . 
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Figure 3: Power law fits for the measured entropy (blue dots) of the strongly interacting Fermi 
gas at 840 G versus its total energy, showing a transition in behavior. Red solid lines show the 
fitted power laws below and above E c — E = 0.41 Ep. Dotted black lines show the extended 
fits. Note that the fit function does not model the smooth transition in slope near the critical 
energy E c , as required for continuity of the temperature. 
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